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Abstract 
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°^ . 

^\ I A spatially discrete version of the general kink-bearing nonlinear Klein- Gordon model in (1 + 1) di- 

mensions is constructed which preserves the topological lower bound on kink energy. It is proved that, 
' provided the lattice spacing h is sufficiently small, there exist static kink solutions attaining this lower 

D , bound centred anywhere relative to the spatial lattice. Hence there is no Peierls-Nabarro barrier impeding 

the propagation of kinks in this discrete system. An upper bound on h is derived and given a physical 
interpretation in terms of the radiation of the system. The construction, which works most naturally when 
the nonlinear Klein-Gordon model has a squared polynomial interaction potential, is applied to a recently 
proposed continuum model of polymer twistons. Numerical simulations are presented which demonstrate 
that kink pinning is eliminated, and radiative kink deceleration greatly reduced in comparison with the 
conventional discrete system. So even on a very coarse lattice, kinks behave much as they do in the contin- 
' uum. It is argued, therefore, that the construction provides a natural means of numerically simulating kink 

dynamics in nonlinear Klein-Gordon models of this type. The construction is compared with the inverse 
method of Flach, Zolotaryuk and Kladko. Using the latter method, alternative spatial discretizations of 
the twiston and sine-Gordon models are obtained which are also free of the Peierls-Nabarro barrier. 

00 

1 Introduction 

A major difficulty in the study of soliton dynamics in the context of high energy physics is that the require- 
ment of Lorentz invariance appears incompatible with integrability, at least for Lagrangian field theories in 
spacetimes of realistic dimension. The PDEs of most interest, therefore, are nonintegrable, and there is no 
prospect of finding interesting exact solutions of multisoliton initial value problems. This is why numerical 
simulation is a popular and important tool. In performing numerical simulations one is forced to discretize 
space in some way, and this inevitably introduces fictitious discretization effects, which one should seek to 
■ minimize. 

The standard discretization of a field equation (in which spatial partial derivatives are replaced by simple 
difference operators, nothing else being changed) replaces it with an infinite system of coupled ODEs, repre- 
senting a network of identical oscillators, nearest neighbours being coupled by springs. Such networks have 
been extensively studied in recent years (l], |3| , primarily because of their condensed matter and biophysical 
applications. The crucial discretization effect encountered is that static solutions may no longer be centred at 
an arbitrary position in space, but must lie exactly on a lattice site, or at the centre of a lattice cell. These 
two types of static solution have different energies (one is a saddle point, the other a minimum), so there is an 
energy barrier, called the Peierls-Nabarro (PN) barrier, resisting the free propagation of solitons form cell to 
cell. As a soliton moves through the lattice, its motion in and out of the PN potential excites small amplitude 
traveling waves in its wake (radiation, or phonons) which drain its kinetic energy, causing it to decelerate. 
Sometimes the kink slows down so much that it has insufficient energy to surmount the PN barrier, whereupon 
it becomes pinned to a lattice cell. 

Clearly, this behaviour is very different from soliton dynamics in the continuum, so if the standard dis- 
cretization is to be used the height of the PN barrier must be made negligible by using a very fine spatial 
lattice, with spacing h of order (soliton width)/20, say. This is computationally expensive. An alternative 
approach, due to Richard Ward, is to exploit the non-uniqueness of the discretization process to find a dis- 
cretization with no PN barrier at all. The hope would then be that the discrete soliton dynamics would closely 
mimic its continuum counterpart even on very coarse lattices, with say h of order soliton width. 
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How does one find such a discretization? Ward's idea is to construct a discrete system which naturally 
preserves the "Bogomol'nyi" properties of its continuum counterpart, that is, a topological lower bound on 
soliton energy which is attained by solutions of some sort of first order self-duality equation. Of course, the 
continuum system one starts with must be of Bogomol'nyi type, but most systems of interest in particle theory 
are (e.g. BPS monopoles, instantons, sigma model lumps, vortices, kinks). The idea has proved very successful 
in the cases of sine- Gordon and </> 4 kink dynamics [Q ^|; that is, it was found that kink pinning was eliminated, 
and radiative deceleration greatly reduced in comparison with the standard discretizations. Unfortunately, 
work on higher dimensional systems (where the benefit of a coarse lattice would be greatest) has been less 
encouraging Q. It is possible to find discretizations of some planar models which preserve the Bogomol'nyi 
bound on soliton energy (and this ensures that the solitons are more stable than they would otherwise be), but 
the bound is unattainable, and there is no discrete self-duality equation. So the PN barrier persists, although 
it is reduced to some extent. Alternatively, by imposing rotational symmetry one can find discretizations with 
saturable Bogomol'nyi bounds but then, of course, the system is effectively one-dimensional. 

The purpose of this paper is to demonstrate that Ward's idea works completely generally for relativistic 
one-dimensional kinks. More precisely, we will construct a "topological" discretization of a general nonlinear 
Klein-Gordon theory in R 1+1 . This will, by construction, preserve the familiar Bogomol'nyi bound on kink 
energy. We will then prove that there exist static kinks saturating this bound centred anywhere relative to 
the lattice, so that there is no PN barrier in this discretization. The construction is a generalization of the 
topological discrete 4 system previously mentioned [^). As an extra example, we apply the construction to a 
recently proposed continuum model of polymer twistons |9) , and again find that elimination of the PN barrier 
leads to continuum-like dynamics deep within the discrete regime. 



2 Relativistic kinks 



The nonlinear Klein-Gordon system consists of a scalar field <f> : R 1+1 — > R whose dynamics is governed by 
the Lagrangian L — Ek — Ep, the kinetic and potential energy functionals being, respectively, 



E K = I 



Ep 



dx i 



dx 



in 



The field equation is the Euler-Lagrange equation obtained by demanding that <f> be a local extremal of 
the action S = J dtL. The following restrictions are made on F c so as to ensure that kink solutions exist: 
F c e C 1 (R,R), 2 < card( J F ( T 1 (0)) < K , and for all u £ F~ l (0), F' c {u) ^ 0. So the model has at least 2 and 
at most countably many degenerate zero vacua, which we label so that F~ 1 (Q) = {ui : i € J} where J C Z is 



a (perhaps unbounded) set of consecutive integers, £ J, and i > j 



> Uj. Given these restrictions, F c 



must change sign at every Uj. It is convenient to label the vacua so that F c is positive on (iii, itj+i) whenever 
i is even (hence negative when i is odd). For any pair i,j 6 J, \i — j\ = 1, by a type kink we mean any 
finite energy solution of the field equation with the boundary behaviour 



<j>{x,t) 



Ui x — ► — oo 
Uj x — > oo, 



(2) 



for all t. It would be conventional to call these "antikinks" in the case where i is odd, but the distinction is 
rather cumbersome for our purposes, and will not be made. 

Static type kinks may be obtained by a Bogomol'nyi argument JTol . To this end, let 



G{cj>) 



d^F c (<j>). 



(3) 



Clearly, G £ C 2 (R, R) and G' = F c . Stable static solutions are local minimals of Ep[<j)\. Assuming 
has boundary behaviour (0), we have, 



E P {4>] 



dx 



dx 



dx^FM 
dx 



2 



> (-1)' / #F C (0) 

J Ui 

= \G( Uj )-G( Ui )\. (4) 

The factor (— l) 4 is included to ensure an optimal bound whether i is even or odd. So, assuming type (i, j) 
boundary behaviour, Ep > \G(uj) — G{ui)\, with equality if and only if the Bogomol'nyi equation, 

t = { - l)lF ^ (5) 
holds. This first order ODE is readily reduced to quadratures. Let Zy be the open interval between Ui and 
Uj. Then, for each (f>o £ ly, there exists a unique solution of (|J) with boundary behaviour (||) and 0(0) = 4>q. 
The set of such solutions constitutes a translation orbit of any one of them, since for each fixed x, 4>( x ) is a 
monotonic function of 4>$. So we may parametrize the moduli space of static kinks, M?j, by 0o £ iij, 
just as well as by kink position, defined for example by b = 4>^ 1 {^{ui + Uj)). The former parametrization 
generalizes more naturally to discrete systems than does the latter. 

For purposes of comparison with the discrete case, we remark that the full field equation is 

4>=^-F^)F c (<p). (6) 
Note that the static field equation is a second order ODE, while the Bogomol'nyi equation (||) is first order. 



3 Topological discretization 



From now on, x takes values in a lattice of spacing h > 0, so x £ hZ = {0, ±h, ±2h, . . .}. We introduce 
the notation /+ and /_ for forward and backward shifted versions of any function / : hZ — > R, that is 
f±(x) := f(x ± h). Also, we denote by A the forward difference operator: Af := h~ 1 (f + — /). Discretization 
of the nonlinear Klein-Gordon system proceeds by defining discrete versions of Ek and Ep. These take the 
form 

h 



Ek 



xEhZ 

where F is a function of <p, </)+ which will be defined below. This will have the correct continuum limit provided 
that lim^^ F{<f>, (/)+) = F c ((j>). 
We shall choose F as follows: 

GQM-G(0) 

1 2 & t 9+ 

Flo.., + ) = { (8) 

F c {<t>) = 0+. 

This looks singular when <f> — 0+j but in fact F : R 2 — > R is C 1 as may be shown by using local compactness 
of R and Taylor's Theorem. (The practicalities of evaluating F in a computer algorithm will be discussed in 
section B.) The important point is that for all 0, <f>+, 

(A<t>)F = AG, (9) 

and this allows one to make a Bogomol'nyi argument analogous to that in section |[ Namely, assuming that 
4> has type boundary behaviour, 

Ep[4>] = ^J2[Aj>-(-iyF] 2 + (-iyhJ2mF 

> (-1)^^ AG 

= \G( Uj ) - G( Ui )\, (10) 
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so this gives a lower bound on type kink energy, which is attained if and only if 



A0 = (-1)\F, (11) 

which will henceforth be called the discrete Bogomol'nyi equation (DBE). 

What this argument shows is that if solutions of the DBE exist with the correct boundary behaviour, 
then they arc minimals of Ep within their class, and hence static solutions of the model. Existence of such 
solutions centred anywhere relative to the lattice will be proved in section ^. Since all such solutions have 
energy \G(uj) — G(ui)\, the kinks experience no PN barrier. 

The discrete field equation is again the Eulcr-Lagrange equation for the action S = / dt (Ek — Ep): 

~ = + -2^ + ^_ _ m 0+) _ pfa (12) 

where Fi is the partial derivative of F : K x M — > M with respect to the first factor, and symmetry of F (under 
cj> «-» (j>+) has been used. The static field equation is a second order nonlinear difference equation, in contrast 
with the DBE ( |i"l| ) which is first order. 

The system should be compared with the standard discrete nonlinear Klein-Gordon system whose kinetic 
energy is the same, but whose potential energy is 

£ P = ^[(A0) 2 +F c (0) 2 ]. (13) 

So the topological discretization differs only in that its "substrate" potential has been symmetrically distributed 
over pairs of nearest neighbours. In the standard discrete system, the Bogomol'nyi argument is lost. The 
discrete field equation is 

^=^-W). (14) 



4 Continuity of the kink moduli space 

It is convenient to make the following definition: a type static kink (where \i — j\ = 1, as before) is a 
monotonic mapping <fi : KL — > K having boundary behaviour (||) and satisfying the DBE ( |ll| ) everywhere. 
Since the continuous translation symmetry of the continuum model has been broken to symmetry under integer 
translations by h, there is no reason to expect a priori, the existence of a continuous "translation orbit" of 
static kink solutions. Indeed, in conventional discretizations, the kink moduli space (that is, the set of 
type static kink solutions) is Mij(h) = |Z, the PN barrier being the difference in E P between "odd" and 
"even" members of this set. By contrast, we will prove that the topological discrete system defined in section 
|]has static kinks centred anywhere relative to the lattice, so M i3 (fe) = R, at least for h sufficiently small. 

The idea behind the proof, which will be presented only for the case of even i (the odd i case follows 
mutatis mutandis), is that any pair (cj),(f> + ) satisfying the DBE ( |Tl| ) may be associated with a point on the 
intersection curve of the plane cj> — 4>+ + hz = with the surface z — F{4>, 4>+) in (0, <f> + , z) space. If h = 0, 
the plane is vertical, and the vertical projection of this curve onto the (0, </>+) plane is clearly the graph of a 
function g interpolating between (m, Ui) and (uj,Uj) (namely, g = Id). If h > 0, but very small, the plane tilts 
slightly away from the vertical, but one may show that such a function persists, and the sequence generated by 
iteration of g on any O S Fj, a solution of the DBE by construction, has the required convergence properties. 



Theorem For any pair i,j £ J such that \i — j\ = 1, there exists /i* > such that for all h £ (0, h*) and for 
all <po £ Fj there exists a type static kink with 0(0) = (fto. 



Proof: Let i be even. Given any C 1 function m : R n — > K, m p denotes its partial derivative with respect to 
the p-th factor. Also, any bounded open interval (z — e, z + e) will be denoted B e (z). We first note that the 
only fixed points of the DBE 

^ = (15) 
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are the vacua, since <f> + = (p =>■ F(cj), <p) — F c ((f)) — 0. Now, one may regard a pair (</>, </>+) satisfying ( |15| ) as a 
zero of the C 1 function / : R 3 -> R 



f{<t> + ,<j>,h) = ct> + -<t>-hF{<t>,ct> + ). (16) 

Note that for any fel, V'j 0) — an d /l(^) "0> 0) = 1 7^ 0, so by the Implicit Function Theorem, for every 
ip G R there exists e,/, > and a C 1 function 17 : B 6i> (ip) x i? eVj (0) — > R such that + = g^,((f>, ft) is the unique 
solution of /(</>+, (f),h) = in B eii (tp) for any (</>, ft,) G domg,/,. Consider any closed bounded interval I Ci. 
Since I is compact, the open covering {B tip {^) : ip 6 7} has a finite subcovering, {B en {if) n ) : n = 1, 2, . . . , A^}, 
which can be used to construct a C 1 function j:/xB f (0)^K where e = inf„ e„ > 0. Namely, 

g(<l>, h) = g^ n ((f>,h) if (f>e(ip n -e n ,ip n + e n ). (17) 

This is well defined by the local uniqueness of each g^ n . For each (cf>,h) € I X (— e,e), + = g(<j),h) is the 
locally unique solution of (|l5|). 

Consider the case I = lij, the closure of the interval Jy. We claim that there exists ft, G (0, e) such that for 
every ft G [0, ft*), g(-, ft) : iy — > R is strictly increasing. To see this, note that g(-, 0) = Id, so <?i(0, 0) = 1 > 0. 
Now 51 is uniformly continuous on, say, 7y x [— |, |], so there exists ft, G (0, |] such that for all ft G [0, ft*) 
<7i(-, ft) : iy — > M is bounded away from 0. 

Since the endpoints Ui, Uj are fixed points of g, an immediate consequence is that for all ft G (0, ft*), g(-, h) 
is an increasing C 1 mapping 1^ — > , which we denote by simply g for the rest of this proof. Hence, for 
any 0o G lij the double-sided iteration sequence 4>(nh) = g n (4>o) exists and is a solution of the DBE ( |l5| ) 
(g° := Id and for n < 0, g n := (g -1 )'"', the existence of the inverse being guaranteed since g is increasing). 
By monotonicity of g, for fixed n G Z, <j>{nh) is an increasing function of 4>q. There can be no fixed points of 
g other than the endpoints (because \i — j\ = 1), so either 

case A: g((f)) > <fi V</> G Zy 
or case B: g((f)) < cj) V0 G lij. 

Given our definition of F c , case A occurs when i < j and case B when i > j (recall we have assumed that i is 
even) . 

Case A: the iteration sequence is increasing, bounded above by Uj and below by Ui and hence both left 
and right convergent. Since Jy is closed, it contains L,R, the left and right limits. Consider the image of 
the forward sequence 4>(nh), n > 1 under g. Since g is continuous, g(4>(nh)) — > g(R). But g((f>(nh)) is a 
subsequence of 4>{nh) and hence converges to R. Hence R is a fixed point of g, and so R — Uj. Similarly, 
considering the image of the sequence 4>(—nh), n > 1, under g _1 one concludes that L — u.i. Hence cj) has the 
correct limiting behaviour. 

Case B: now the iteration sequence is decreasing and bounded above by Ui and below by Uj. A similar 
argument shows that (j)(nh) — > Uj as n — » — 00 and <f>(nh) — > itj as n — > 00. This completes the proof. □ 

An important point to note is that ft*, the upper bound on lattice spacing up to which My (ft) is continuous 
(henceforth assumed to be optimal), depends in general on (i, j). If i r c ~ 1 (0) is finite, a global upper bound on 
ft may be found, namely, 

h*:=wf{K(i,j):i,je J,\i-j\ = l}>0 (18) 

so that the topological discretization with < ft < ft* has continuous moduli spaces for all types of kink. This 
is not necessarily true if i r c _1 (0) is infinite (for then it is possible that ft*" = 0). For any particular one 
can find an upper bound on ft* however, meaning that in no case does continuity of My (ft) persist for all 
ft > 0. 

The upper bound is obtained by again thinking of the graph of the function g(-,h), whose existence was 
proved above, as the vertical projection of the intersection of the plane 4> — <f>+ + (— l) l hz — with the surface 
z = F((j) 1 4>+)- Recall that this graph always passes through (uk,Uk), where k = i,j. Consider the behaviour 
of the tangent to the graph at (iik, Uk)- As ft grows large, this tangent passes either through the horizontal or 
through the vertical (depending on whether i is odd or even, and whether k — i or k = j). In either case, from 
this point on (i.e. for larger ft) either monotonicity or existence of <?(•, ft) is lost, and so the existence or correct 
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limiting behaviour of the iteration sequence 4>(nh) = g n (cf)o,h) is generically lost. This will be explained in 
more detail shortly. 

The gradient of the tangent through {u k , Uk) may be calculated by demanding that (</>, <f> + ) = (u k + S,Uk + 
S+) remains a solution of the DBE ( |ll| ) to leading order in ||<S|| = [S 2 + Using Taylor's Theorem for G, 



one finds that 



G(u k + S) = G{u k ) + -F' c {u k ) + o(5 2 ), 



-lYhF^u k ) 



i-iyhF' c (u k ) 



(19) 



(20) 



So if i is even, meaning that F' c (ui) > and F' c (uj) < 0, the tangent through (ui,Ui) is vertical when 
h = 2/F' c (ui), and the tangent through (uj,Uj) is horizontal when h = 2/\F' c (uj)\. On the other hand, if i is 
odd, then F' c {ui) < and F' c (uj) > 0, so the tangent through (ui,Ui) is horizontal when h = 2/|F c '(ui)|, and 
the tangent through (uj,uj) is vertical when h = 2/F' c (uj). In either case, we obtain the same upper bound 
on /i*, namely 

<hi ~^ItJ-^,ttJ-tj)' (21) 



\F>( Ui )\> \Fl( Uj )\ 



One may similarly define 



hi := sup 



(22) 



\Fl(ui)\ '|^(%)l 
which will prove useful below. 

A typical situation is depicted graphically in figure 1. Here, hi = 2/|F c '(ui)|, h 2 — 2/|i^(u )| and hi < h < 
/i2, so the tangent through (ui,ui) has passed through the horizontal. It is impossible for cj)(nh) to converge 
monotonically to ui with ui as a cluster point of the sequence: the sequence must overshoot ui as 4>(nh) 
exceeds tp h . The only possibility for a (0, 1) kink to exist for this value of h, therefore, is that </> has a constant 
right hand tail, that is, there exists TV € Z such that (j>(nh) = ui for all n > N. This, in turn, is possible only 
if 0(0) = ifh, or one of its preimages under the iteration. One can restate this condition as (f>(0) £ where 



*h== {ih(.<Ph):n = 0,1,2,...}, 



(23) 



7ft, being the partial inverse of <?(-, h), with 7^ : [uq, u{\ — > [uo, V/i]- The set $ft parametrizes a discrete moduli 
space of kinks Moi(/i), all identical modulo lattice translations. Note that the left hand boundary behaviour 
of <p remains good for all 4>(0) G (mq,mi) in this situation. 




Figure 1: A section of the curve A(j) — F, with hi < h < h 2 . 



In figure 2, h has been increased further, so that h > h 2 , and now the sequence 4> behaves badly at both 
ends. In order for (j)(nh) — > uq monotonically as n — » —00, (f> must have a constant left tail, so one must have 
0(0) <E $ft, where 

$ fc := {TOO :n = 0,1,2,...}, ( 24 ) 
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and "fh ■ [uo,tph] ~ y VPh^i] is denned by the curve shown. In this case, for (j) to satisfy the definition of a 
static kink, both the left and right tails must be constant, which can only happen if $^ = One expects 
this to occur only for a discrete set of h values, which accumulate towards h 2 . When h becomes very large, 
one expects the projected intersection curve to pass outside the open square (tto, u\) x (uq, u\) completely, so 
that no static kink solutions are possible. 




Figure 2: A section of the curve Aip = F , with h > hi 



We are led, therefore, to the following conjecture for the generic behaviour of Mij(h) as h varies: for 
h E (0, hi], Mij(h) is continuous, M; for h 6 [hi,h 2 ], Mij(h) is discrete, KL\ for almost all h 6 (h 2 ,oo) 
Mij(h) — 0, but there exists a bounded countable set H C (h 2 , oo) with inf H — h 2 such that My(/i) = hTL 
for all h £ H. The conjecture is represented schematically in figure 3. It matches the observed behaviour 
for the topological discrete sine-Gordon j| and 4 || systems. In both these cases hi = h 2 because the 
kinks interpolate between identical vacua, so the middle band is missing, and for the sine-Gordon model 
H = 0. The example of twistons, discussed in section produces the full generic behaviour, with all three 
bands. It is probably possible to construct a function F c such that this conjectured behaviour fails (by making 
F c {\{ui + Uj)) very close to zero, for example), but such a system is rather contrived. One should note that 
the onset of discreteness in the kink moduli space is due to the requirement that kinks should be monotonic 
(after all, they are supposed to model continuum kinks, which always are), rather than the appearance of 
some kind of PN barrier. 

h 





h=h 1 



h=h. 



I (kink position)/h 

Figure 3: Schematic picture of the generic structure of the kink moduli space Mij (h) as h varies. Horizontal 
sections represent the moduli spaces for each specific value of h. 

The upper bound /i» < hi has a remarkable physical interpretation in terms of the radiation of the system. 
Just as different pairs of vacua Ui, Uj have different types of kink interpolating between them, so each different 
vacuum Uk has its own type of radiation. In each case, the dispersion relation for this radiation may be 
obtained by substituting a traveling wave ansatz, cos(kx ~ Lut) , into the discrete field equation linearized about 
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4> = itfe. Noting that F\{uk, Uk) = ^F' c {u^), one finds that u> and k must be related by 



= [K{u k )f 



~ 2 {K{n k )f 



sin 



kh 

T 



(25) 



Note that this sinusoidal dispersion relation collapses to a flat line when h = 2/\F' c {uk)\ 1 and that for larger h, 
the curve has a maximum instead of a minimum at k = 0. For h > 2 / li 7 ^ ) | , then, the discrete system fails 
to model accurately the radiation of the continuum system (which obeys the relativistic energy-momentum 
relation ui 2 = [F' c {uk)\ 2 + k 2 ) even in the long wavelength limit. So the upper bound on h* derived above 
is precisely that value of h for which radiation about one or other of the vacua between which the kinks 
interpolate starts to behave badly. 



5 Twistons 

The original motivation for the construction in section || was to provide an efficient and natural means of 
simulating soliton dynamics numerically on a computer. In this light, two substantial objections to the 
discretization procedure can be raised. First, even if F c is known in closed form, its anti-derivative G may not 
be. In this case, every evaluation of F in the computer program would require the approximate evaluation 
of two definite integrals (G((f>) and G((f>+)). Clearly this cannot be computationally efficient. Second, the 
piecewise definition of F is highly inconvenient for numerical purposes. That is, even if an explicit formula for 
G exists, evaluation of F for close to </>+ would be subject to significant rounding errors, unless (for example) 
a polynomial approximation to F were used in a narrow strip containing the diagonal (j) = </>+. Again, this 
is complicated and inefficient. There is, however, a large class of potentials F 2 for which both objections are 
avoided, namely those where F c is polynomial. In this case, G is also polynomial (hence available in closed 
form) and, by the remainder theorem, so is [G((j> + ) — G(0)]/(0+ — <j)), which may obviously be extended to 
the whole (</>, 0+) plane to give a global, rather than piecewise, definition of F. 

One example of such a polynomial discrete system has already been investigated ||. To illustrate the 
construction further, we shall investigate the topological discretization of a recently proposed continuum 
model of so-called twistons in crystalline polyethylene || . We emphasise that the aim here is not to propose a 
physically realistic discrete model of such polymer twistons, but rather to demonstrate that our discrete system 
performs favourably in comparison with the conventional discretization when simulating soliton dynamics in 
the continuum model. 

The continuum system has interaction potential F 2 , where 

FM =0(l-0 2 ). (26) 

We have normalized F c differently from and absorbed a coupling constant. So, there are three zero vacua, 
Uo = — 1, Ui = and u% = 1, and four different types of kink. Actually, all four types are trivially related to 
one another by reflexions in space and/or the codomain, and so we shall consider only type (0, 1) kinks. The 
static profiles of these may be found explicitly using the Bogomol'nyi argument outlined in section 0, namely, 



i(l-tanh(a;-6)). (27) 

Turning now to the topological discretization of this model (so now x 6 hZ), we find that G(4>) = <f> 2 (<fi 2 — 
2)/4, and 

m 4> + ) = G{ % ] ~ G ^ ] - + <t> 2 2)(0 + + 0). (28) 

As remarked above, F has a simple closed form, and its evaluation presents no problems for a computer 
program. By the general existence theorem of section |[ there exists > such that for all h G (0, /i*) there 
is a continuum of (0, 1) kink solutions Moi(/i) of the DBE, 

^ ± T ± = \{<t>\ + tf-m + +^ (29) 
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parametrized by cj>(0) £ (—1,0). All such solutions have energy Ep — \G(Q) — G(— 1)| = i, so there is no PN 
barrier. Recall that, in general, one has an upper bound hi on h*. In the present case, F' c [— 1) = —2 and 
F' C (Q) = 1, so hi = 1 and ft 2 = 2. 

We do not propose to prove rigorously that h* = hi = 1 for this model. Rather, by plotting the zero 
contours of the function 

/(</>+>, ft) = 0+ - <j> - + 2 - 2)(0+ + 0) (30) 

for various values of ft. we will demonstrate that the generic behaviour of My (/i) described in section ^ is highly 
plausible. Recall, from the proof in that section that this contour contains, at least for small h, the graph of the 
function g(-, ft) which one iterates to generate the static kink. Zero contours of / for ft, = 0.5, 1.0, 1.5, 2.0, 2.5 
and 4, found using the contourplot utility of Maple, are shown in figure 4. While ft < 1, the contour contains 
the graph of a continuous increasing function g(-, ft) : [—1, 0] — > [—1, 0], so that Moi(ft) is continuous. However, 
when ft = 1, the tangent to the contour at (—1, —1) passes through the vertical, so for ft > 1 only kinks with a 
constant left tail can exist, which can only happen for nongeneric 4>(0), and M i3 (fe) is discrete. When ft = 2, 
the tangent to the contour at (0, 0) passes through the horizontal. For ft > 2 kinks must also have a constant 
right hand tail, so they may only exist for nongeneric (/>(0) and ft. This situation continues until ft > 4, when 
the contour passes outside the open square (—1,0) x (—1,0). Indeed, at ft = 4, the DBE supports only "step 
function" kinks, 

*M-{~o.">$ (31) 

as may be shown by explicit calculation. 




Figure ^: Zero contours of (f> + — (f) — hF(<f>, 4> + ) for the values ft = 0.5,1.0,1.5,2.0,2.5 and 4, from least to 
most curved. 

Since the DBE (|2^) is a cubic polynomial equation, it is possible in principle to solve for (f> + in terms of 
<j) (or vice versa) and hence generate exact static kink solutions, although it appears that no closed formula 
for 4>(x) exists (in particular, a continuum- inspired ansatz does not seem to work). In general, if F c is a 
degree p polynomial, the DBE is a degree p polynomial equation for + (0), so for p > 5 such an approach 
would not work. Even for p = 3, 4, implementation would be complicated. Instead, one may solve the DBE 
approximately with some kind of root-finding algorithm. In figure 5 we present a kink solution found by 
solving (^9|) with ft = 0.8 using the bisection method. Note that the kink structure is spread over very few 
lattice sites, so the system is deep in the discrete regime. 
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Figure 5: Example of a static kink solution on a lattice of spacing h — 0.8 obtained by approximate solution 
of the DBE using the bisection method. Plot (a) shows the field <j>. Plot (b) shows the energy distribution of 
the kink. The energy should be thought of as located in the links between pairs of lattice sites. This plot shows 
the energy in each link as a fraction of the total energy. 



6 Numerical simulations 

So far, we have considered only the static model. In this section we will discuss numerical simulations of 
the full, time-dependent discrete field equation, that is, numerical solutions of the coupled set of ODEs (|l2|). 
In solving these ODEs one must also discretize time in some way. It should be regarded as a strength of 
our approach that there is great freedom in precisely how one does this, since one may choose a method 
specifically adapted to the dynamics one seeks to simulate. For example, in kink-antikink interactions, the 
solitons usually spend a long time well separated, then experience a brief, but very violent interaction when 
they finally become close. In this situation a variable time step (rather than a fixed, rectangular space-time 
mesh) is clearly sensible. For other problems it may be desirable to preserve the symplectomorphicity of the 
phase space flow (both (JsJ) and (|l2] ) are Hamiltonian dynamical systems) , so that a symplectic integrator with 
fixed time step would be more appropriate. 

Since the aim here is to test our discretization's performance, we shall consider a simple and clear-cut 
dynamical problem, for which the continuum behaviour is easily understood: single (0, 1) kink motion at 
constant speed. The continuum field equation (0) is Lorcntz invariant, so there exist Lorentz boosted kinks 
moving with any speed v € [0, 1). These provide the initial data, sampled on the discrete lattice hZ, for our 
simulations. The lattice spacing was chosen to be h = 0.8, so that the kink width is comparable to h. Recall 
that a major departure of conventional discrete systems from their continuum counterparts is that solitons 
suffer radiative deceleration, and possibly pinning, in the highly discrete regime. An important test of the 
topological discretization, therefore, is to monitor kink velocity v(t) over simulations of long duration (3000 
time units in this case), to see to what extent radiative deceleration is present. One certainly expects moving 
kinks to excite some radiation. For a type kink, this will move with a maximum group velocity of 

= l-~l^(«*)l (32) 

where k = i if v > or k = j if v < 0, since phonons are excited behind the kink, and hence arc small 
oscillations about its trailing vacuum |fj"l|| . So for right moving (0, 1) kinks, «™ ax = 0.2, while for left moving 
(0, 1) kinks, i;™ ax = 0.6. In either case, one wishes to avoid radiation being reflected from the fixed boundaries 
and interfering with the kink's motion, so the first few sites at either end of the (finite) lattice are damped. 
This adds a non-Hamiltonian piece to the flow, so a symplectic integrator is not appropriate, and we choose to 
use a fourth order Runge-Kutta method with fixed time step St = 0.01. The algorithm, and the checks made 
on its accuracy and stability (energy conservation etc.) have been described previously ||. 

Simulations were performed for various initial velocities, and in all cases it was found that kink pinning was 
eliminated, as one would expect: kinks propagate freely, indefinitely. Figure 6 shows plots of v(t) for right- 
moving kinks with initial velocities v(Q) — 0.2, 0.4, 0.6 and 0.8. Figure 7 presents similar graphs for left moving 
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kinks, with v(0) = —0.2, —0.4, —0.6 and —0.8. In every case, the upper curve shows v(t) for the topological 
discretization (112), while the lower curve shows v(t) for the same initial value problem for the conventional 
discretization (fL4), solved using the same numerical scheme. The thickness of the curves is due to velocity 
oscillations as the kink moves from cell to cell through the lattice. These are partly an artifact of the way 
we define kink position (linear interpolation in this case), and partly dynamical. In both systems one may 
regard the effective kink mass as depending periodically on the kink position, and in the conventional (but not 
the topological) discretization, the kink velocity oscillates as it moves up and down in the PN potential. As 
found previously for sine-Gordon and 4 kinks, radiative deceleration is drastically cut using the topological 
discretization. 
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Figure 6: Comparison of the motion of right-moving kinks in the topological discrete twiston system with that 
in the conventional discretization, for various initial velocities. The lattice spacing is h — 0.8. In each case, 
velocity is plotted against time for both systems on the same graph, the lower curve showing the data for the 
conventional discrete system. 
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Figure 6: Comparison of the motion of left-moving kinks in the topological and conventional discrete twiston 
systems, with lattice spacing h = 0.8. Speed is plotted against time for both systems on the same graph, 
the lower curve showing the data for the conventional discrete system. As with right-moving kinks, radiative 
deceleration is much less pronounced in the topological discretization. 



7 Other methods 

The discretization method outlined in section [| is not the only way to arrive at discrete systems without a 
Peierls-Nabarro barrier. Having performed the construction, one is always free to replace A<fi by D = XA(f> 
and F by F — F/X where X is any function of (4>,(f)+,h) which reduces to unity in the continuum limit, 
h — » 0, <fi+ = 4> + O(h). Since DF = AG, the potential energy functional 



E, 



(33) 



is again subject to the the BogomoPnyi argument, yielding the same lower bound on type kink energy, 
but the DBE is now different: D = ±F. For example, the sine-Gordon model has F c — sin, and hence its 
topological discretization has 



n<t>, 



cos <p + — COS ( 

<j>+-(j> 

sin d> 



7^ 



(34) 



If one chooses X = sa\c^{<j) + — 0), (where sincx := x 1 since when i / and sincO := 1), then defining 
D = XAcj) and F — F/X one recovers the original topological discrete sine-Gordon system [||, namely 



D 



F 



— sin — 

h 2 



(35) 



In general, the kind of analysis carried out in section |^ should still work for the new DBE (and certainly 
does in the specific case of the sine-Gordon system outlined above), although the geometric picture of a plane 
intersecting a surface no longer applies. Using this trick, one may be able to avoid a piecewise definition of F 
even when F c is not polynomial. However, no systematic method for choosing X suggests itself. 
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Finally, we wish to describe briefly a completely different approach due to Flach, Kladko and Zolotaryuk, 
called the inverse method. The following is a somewhat reinterpreted version of the method outlined in fl^] 
(the authors of which were primarily concerned with exact moving discrete solitons, rather than the Peierls- 
Nabarro barrier). The idea is that, given a continuum nonlinear Klein-Gordon model possessing a type 
kink with static profile <j>x, one seeks a one-site substrate potential Vh(4>) such that the discrete field equation 

4> = ^ V h (4>) (36) 

preserves the one parameter family of static solutions <p(nh) = <pK{nh — b). Clearly, one needs 

V^k(z)) = ^M>k{z + h)- 2<t> K {z) + 4>k{z - h)} (37) 
hr 

for all zel. This uniquely determines Vl : l t j —* R by monotonicity of 4>k- Note that, by construction, in 
the limit h — > 0, the right hand side of (B7j) becomes <$>" K so, given that <\>k satisfies the static continuum field 
equation (^[), V^(</>) must have the correct continuum limit (lim/^o = F C F' C ). Since this is a continuous 
family of static solutions, they must all have the same energy. Indeed, equation ([?7|) may be interpreted as 
the condition 

d -E P [4> K {nh-b)] = 0, (38) 



db 

where the potential energy is 

E P [4>] = hJ2 



(39) 



Hence there is no PN barrier resisting the propagation of kinks in this discrete system. If tpK is analytic, 
Vl((f>) can be written as a power series in <f>, although there may exist no closed formula for V£. This may be a 
significant drawback of the inverse method, in comparison with the topological discretization approach, which 
works very explicitly, at least for all polynomial F c . 

In some special cases, can be obtained in closed form: one example is <fr theory, where 4>k = tanh 
and an explicit formula for Vh has been known for some time pljfl . As a new example, one could consider the 
continuum twiston model defined in section ^, whose (2, 1) kink has profile 4>k{%) = [5(1 — tanh x)] 2 . In this 
case, 4>k{z ± h) can be written in terms of 4>k(z) using hyperbolic trigonometric identities, so one finds, 



... <t> J / 1 — tanh/i / 1 + tanh ft, I 

VhW = 7? jy 1 + (1 - 202)tanhft + V 1 - (1 - 2^)tanh/ l " 2 J ' (40) 

which has W(±l) = Vu(0) — 0, as expected. It is unclear how "special" F c must be in order for the Vl obtained 
by this inverse method to be expressible in closed form - perhaps this happens whenever F c is polynomial. 
If so, one would expect the expressions involved to become extremely complicated as the degree of F c grows. 
The inverse method also works explicitly when applied to the sine-Gordon kink profile, 4>k{x) — 2 tan -1 e x , 
so generating yet another discrete sine-Gordon system with no PN barrier (see also flij ). In this case, 



e h tan - I + tan -1 ( er h tan - 
2 J \ 2 



(41) 



Equations ( p0[) and (^Tj) can be integrated to give explicit formulae for Vh, and hence Ep, for these models, 
although the formulae are rather complicated. This is little practical disadvantage, since only appears in 
the discrete field equation. The strength of this method is that one has an explicit formula for the family of 
static discrete kinks, instead of an abstract existence theorem. It certainly merits further investigation. 



Acknowledgments 



The author wishes to thank S. Flach and Y. Zolotaryuk for their patient explanation of the inverse method. 



13 



References 

[1] M. Peyrard and M.D. Kruskal, "Kink dynamics in the highly discrete sine-Gordon system" Physica D14 
(1984) 88. 

R. Boesch, C.R. Willis and M. El-Batanouny, "Spontaneous emission of radiation from a discrete sine- 
Gordon kink" Phys. Rev. B40 (1989) 2284. 

D.B. Duncan, J.C. Eilbeck, H. Feddersen, J. A. Wattis, "Solitons on lattices" Physica D68 (1993) 1. 

J.M. Speight and R.S. Ward, "Kink dynamics in a novel discrete sine- Gordon system" Nonlinearity 7 
(1994) 475. 

J.M. Speight, "A discrete (j) 4 system without a Peierls-Nabarro barrier" Nonlinearity 10 (1997) 1615. 

R.S. Ward, "Bogomol'nyi bounds for two-dimensional lattice systems" Commun. Math. Phys. 184 (1997) 
397. 

R. Leese, "Discrete Bogomol'nyi equations for the nonlinear 0(3) sigma model in (2 + 1) dimensions" 
Phys. Rev. D40 (1989) 2004. 

T. Ioannidou, "Soliton dynamics in a novel discrete 0(3) sigma model in (2 + 1) dimensions" Nonlinearity 
10 (1997) 1357. 

D. Bazeia and E. Ventura, "Topological twistons in crystalline polyethylene" preprint cond-mat/9810249"[ 

E. B. Bogomol'nyi, "The stability of classical solutions" J. Nucl. Phys. 24 (1976) 449. 
J.M. Speight, Ph.D. Thesis, Durham University (1995), pp 47-48. 

S. Flach, Y. Zolotaryuk and K. Kladko, "Moving lattice kinks and pulses: an inverse method" preprint 
1998. 

V.H. Schmidt, "Exact solution in the discrete case for solitons propagating in a chain of harmonically 
coupled particles lying in double-minimum potential wells" Phys. Rev. B20 (1979) 4397. 

W.J. Zakrzewski, "A modified discrete sine-Gordon system" Nonlinearity 8 (1995) 517. 



14 



